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The variational principle of V. I. Arnold [J. Appl. Math. 
Mech. Vol. 29, P. 1002 (1965)] is extended to the general con- 
servative inhomogeneous, compressible, and conducting fluid. 
The concept of iso-vortical flows is generalized to an "invari- 
ant foliation" of the phase-space. The foliation, which may or 
may not correspond to explicit conservation laws, is derived 
from the equations of motion and used for Lyapunov stabil- 
ity. A nonlinear three-dimensional (magneto-) hydrodynamic 
stability criterion is formulated. 



The standard approach to hydrodynamic stability in- 
volves linearization about an equilibrium flow in order 
to solve for eigenfrequencies or establish a Lya- 

punov stability criterion for the linearized system 
A variety of linear variational principles were developed 
for both neutral fluids and magnetohydrodynamics 
(MHD) in which the stability criter ion is expressed 

in terms of a positive definite quadratic form. It is well 
known that the linearized stability does not guarantee 
the true (Lyapunov) stability, such as in the toy system 
du/dt = u^, whose equilibrium u = is linearly stable. 

Nonlinear stability is guaranteed by the presence of 
an integral of motion, for example, the energy H, which 
assumes a non-degenerate extremum (a minimum or a 
maximum) subject to the conservation of any other inte- 
grals of motion, for example, Casimir invariants The 
possibility to write explicitly a full infinite set of integrals 
is mostly limited to two-dimensional systems. By explicit 
we mean an integral of motion which can be written in 
terms of the physical flelds of velocity, density, etc., in a 
way which does not require the solution of the equations 
of the motion. In three dimensions, such integrals are 
scarce. For example, the Euler equation 

dtu^ = V X (v X Lj), a; = V X V, V • v = 0, (1) 

conserves explicitly only the energy H and the helicity /: 



= (p -v ■ d£ = / u) ■ dS = const, 



I = 



(2) 



(Here and below, unless specified, all integrals are under- 
stood over the domain occupied by the fluid. An approri- 
ate conservative boundary condition, such as zero normal 
velocity, is implied.) In addition to the explicit integrals 
(0), there is also the infinity of Kelvin invariants. 



(3) 



expressing the velocity circulation around (or the vortic- 
ity flux through) any closed contour 7(f) moving with the 
fluid velocity v. Integrals (|^) are implicit in the sense 
that their definition involves contours 7 whose motion 
must be solved from Eq. (|^). Although there is no appar- 
ent way of incorporating integrals like (^ in a Lyapunov 
functional, Arnold j|] proposed that the conservation of 
all vorticity integrals (H) has the geometrical meaning 
of confining the system to an "iso-vortical sheet" in the 
infinite-dimensional phase space. Different sets of initial 
vorticity integrals specify different sheets such that the 
whole phase space is "foliated," as if by iso-surfaces of an 
integral of motion (Figure 1). 




FIG. 1. The infinite-dimensional phase space of all incom- 
pressible fluid flows, t<j(x), is foliated by iso-vortical invariant 
sheets parameterized by the displacement function ^ (x). Each 
sheet is an inflnite-dimensional subspace of tu. The dynamics 
keeps an orbit on a sheet. 

The usefulness of the foliation for stability is due to the 
local explicit parameterization of the iso-vortical sheets 
by an incompressible "displacement" ^(x), such that vor- 
ticity fields sharing the sheet with the reference flow 
a;o(x) are written uj = u)o + Sujq + ^S^ujq -I- . . ., where 



(5u; = V X (^ X a;). 



(4) 



The iso-vortical variation operator d derives from the 
"modified dynamics" dtu: = V x {dt^ x uj), V • ^ = 0, in 
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which the vorticity is incompressibly advected in a way 
similar to the Euler equation (|l]) , but by the velocity field 
dt^ entirely unrelated to the actual flow v. Since the con- 
servation of the vorticity integrals (3) is independent of 
the relation between v and lj, Eq. (4) follows. 

Arnold variation (^) makes the Hamiltonian H station- 
ary if and only if the flow v is in equilibrium. Then the 
second energy variation, 



S'^H = J ((5v)2 - I X a; • V X (I X v), 



(5) 



if definite for all incompressible ^, guarantees that the 
equilibrium is stable 

Given this long introduction, we briefly report on a 
generalization of the Arnold method in two important 
ways. Firstly, our fluid equations include com- 

pressibility, varying entropy and also magnetic fleld, but 
still no dissipation. In such a general formulation, it is 
difficult to write all integrals generalizing Eq. (||) . There- 
fore, and secondly, an analog of the iso- vortical variation 
is formally derived from the dynamics, without regard to 
either explicit or implicit integrals of motion. A new out- 
come of this procedure is a variational principle for ideal 
MHD stability with fluid flow, a long-standing plasma 
physics problem ||,|,|l]j|]. 

We consider the following hydrodynamic equations for 
an inviscid, ideally conducting fluid: 



p(5tv -I- V ■ Vv) = - Vp(p, s) -I- j X B 
StB = V X (v X B), 
9ts + V ■ Vs = 0, 
9tp + V • pv = 0. 



(6) 
(7) 
(8) 
(9) 



Here p is the fluid pressure, p the density, s the entropy, 
(j) the external gravitational potential, B the magnetic 
field, and j = V x B the electric current. The fluid flow 
conserves the energy 



H = 



P 



Y + <p, 



B2 
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(10) 



where e is the specific internal energy defined by the stan- 
dard thermodynamic relation 



de ^ T ds — pd{l/ p). 



(11) 



The varying entropy and the Lorentz force in Eq. (|^) 
break the "frozen-in" law for the vorticity u; = V x v, 
and intcad of (^ wc now have 



dfiv = V X V X u; 



. B 

-J X — 
P 



P dsP dp . 



(12) 



By analogy with the iso- vortical variation M), one can 
introduce modified dynamics for Eqs. (^)-(^ in many 
different ways |12| . Our choice is dictated by the desire to 



have a zero energy variation for equilibrium fiows. Upon 
a dozen of attempts, the following procedure works to 
our satisfaction: We replace v dt^ in Eqs. (|7|)-(||) and 
(|T^). In (|T^), we also write j dtT] (V • tj = 0) and 
replace the integral by a scalar dta. The result is the 
generalized iso-vortical variation, 

Sv = ^xuj + T]x h aVs + V/3, = V x C, (13) 

P 

(5B = Vx(^xB), Ss^-i-Vs, Sp=-V-p^, (14) 

which depends on two arbitrary vectors ^ and C and two 
arbitrary scalars a and f3. 

Although variation (^^-(^4|) conserves magnetic flux 
and entropy integrals (which are well known and not writ- 
ten here), it is unclear what other conservation laws, if 
any, are accounted for by this variation. Nevertheless, 
the derivation above clearly implies that the phase-space 
sheets parameterized by (^,C)C«,/9) are invariant sheets, 
which can be interpreted as iso-surfaces of some integrals 
of motion and thus used for stability analysis. In the 
limit of zero magnetic field and constant entropy, varia- 
tion dl^ ) reduces to Arnold's iso-vortical variation (^). 

The number of arbitrary functions in the variation 
(p3|)-(p^ by no accident equals the number of dynami- 
cal equations (||)-(||). Upon varying the energy ( p^ ) and 
using Eq. (0), a few integrations by parts yield 



6H= J I • (pv • Vv + Vp - j X B + pV(/)) 

-C • V X (v X B) + apv • Vs - /3V • pv (15) 

By design, the condition that 5H — for all {^jC, 
is equivalent to an equilibrium solution of Eqs. (|3)-(|9|). 
The second variation of velocity, 

S^v ^^xSLj + rjxS h aV5s + V/3, (16) 

P 

and similar expressions for (5^(B,s,p) are now used to 
calculate the second energy variation: 



6^H = S^pe 



2JV 2'*" 

+p{Sw)^ +pw-S^v + 2 6pw-5w. (17) 

Equations (|l3|)-(0 define S^H as a functional of 
(^,Ca,/3). Two comments regarding the form of 6^H 
are in order. 

First, the suspicious linear term V/3 in the second ve- 
locity variation ( p^ ) is multiplied by an incompressible 
pv in Eq. (^) and thus vanishes upon integration by 
parts. So, as it should, S^H is a quadratic functional of 
the independent variables a, (3). 

Second, the integrand of ( p^ can be written as a 
quadratic polynomial of a with the coefficient p(Vs)^ 
in front of . Therefore, the definite sign of 6^H can be 
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only positive, and, for this, it is necessary and sufficient 
that the reduced quadratic form of (|, C, /?), 

-p{n ■ 5"^ff + (5'v • {p5"w + wSp) + pv ■ -q x S—, (18) 

P 

be also positive definite. Here n = Vs/|Vs|, and the 
space-saving notation is introduced, 

(5'v = |xu;+J7xB/p+V/3, S"v = ^'v+^■V^r-v■V$. 



No further "simple" minimization of Eq. (|18|) is possible. 

To make Eq. (|l^) more explicit, the internal energy 
variation can be transformed as 

Js\pe)= Jv-^ipdppV-^ + ^-Vp), (19) 



and the magnetic energy variation 



(20) 



Also, the second density variation is S^p = V • (^V • p^). 

Equation (|l8|) is the main result of this paper. It can be 
used to establish nonlinear stability of general MHD equi- 
libria with fluid flow, for example numerically, by test- 
ing the sign of W for sets of seven scalar trial functions 
(^7 Ci /?)■ Since all known explicit and implicit integrals of 
motion are incorporated in our scheme (the possibly con- 
served linear and angular momenta amount to choosing 
an appropriate frame of reference), we propose that the 
sufficient stability criterion W > is also necessary for 
the true nonlinear stability of an ideal MHD equilibrium. 
This conjecture is also supported by the static limit of 
zero flow, v = 0, in which our variational principle re- 
duces to the sum of Eqs. (|l9| ) and (po|), or the standard 
MHD energy principle whose violation means a lin- 
ear instability. As a by-prodict, we thus find that the 
linear stability criterion of Bernstein et al. ^ for static 
equilibria is also a nonlinear stability criterion. In a gen- 
eral situation with fluid flow, an indefinite W may not 
result in an exponential instability, but rather lead to 
a slower, nonlinear perturbation growth and subsequent 
turbulence. This scenario will be described elsewhere. 

The other two limiting cases we would like to mention 
are (a) the hydrostatic equilibrium with (j) — gz and v = 
B = and (b) the incompressible neutral fluid with p = 
s — const and B = 0. For the former case, the condition 
W > yields the well known convective stability criterion 
ds/dz > and dp/dz < 0. In the Euler limit, the 
incompressibility is introduced by letting the sound speed 
= dpp to infinity. The minimum of Eq. (jl^ ) then 
implies V-^ for the "most dangerous" perturbations, 
and further minimization of (|l7|) with respect to (3 results 



in V • = and the restricted Arnold criterion that 
Eq. (P) be positive definite. 

Thus, in several evident limiting cases, our stability 
criterion reduces to the already known results. How- 
ever, in the general case of flow of an inhomogeneous 
fluid, with or without magnetic field, it appears new. 
The exact mathematical meaning of the generalized iso- 
vortical variation and the status of the resulting stability 
criterion Vl^ > 0, Eq. (p^), remain unclear to this au- 
thor. For example, no a priori estimates exist for three- 
dimensional hydrodynamic perturbations, unlike those in 
two dimensions, where all Casimir integrals are explicit 
|l0|,^. On the "physical level," the sufficient stability 
criterion W > looks rigorous, less so as a necessary 
one. 
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